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Abstract. The extended Gaussian ensemble introduced recently as a generalization of the canonical en- 
semble, which allows to treat energy fluctuations present in the system, is used to analyze the inelasticity 
distributions in high energy multiparticle production processes. 
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1 Introduction 



The high energy multiparticle production processes are 
very important source of information on the dynamics of 
hadronization process, in which some amount of the ini- 
tially available energy is subsequently transformed into 
a number of secondaries of different types. Such processes 
can be described only via phenomenological models, which 
are stressing their different dynamical aspects, like spe- 
cific energy flows £Q or their apparent thermal-like char- 
acter Actually most of the characteristic features of 
hadronization can be described in universal manner by 
means of Information Theory (IT) approach, both in its 
extensive |3| or nonextensive jHEK] versions. The main 
difference between them is that whereas former is using 
only energy- momentum conservation constraint, the later 
accounts also for some intrinsic fluctuations present in the 
hadronization process, either in the form of fluctuations of 
temperature or in the form of fluctuations of the num- 
ber of produced secondaries U 1 . Recently the extended 
gaussian ensemble (EGE) approach has been proposed to 
account for some fluctuations in statistical mechanics and 
it was presented also in the IT formulation |10| . The ques- 
tion, which we would like to address here, is whether EGE 
can find application in deducing some new information 
from hadronic production processes. 
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2 Extended Gaussian ensemble from IT 



Following |3]0]|o]|b] we are interested in applying IT to 
deduce the most probable and least biased energy dis- 
tributions of particles produced in hadronization process 
in which mass M transforms into given number N of 
secondaries of mass \i and mean transverse mass px = 
(pt) 2 each, distributed in the longitudinal phase 



VP 

space described by rapidity variable, y (such that energy 
of particle is E = \xt cosh y). We are therefore interested in 
(normalized) rapidity distribution p(y) = (l/N) ■ dN/dy, 
J dyp(y) = 1, which according to IT |3] is obtained by 
maximizing Shannon entropy 



S 



dyp(y) \np(y), 



(1) 



under condition of reproducing known a priori mean value 
of energy of produced secondaries (K denotes the so called 
inelasticity of reaction to be discussed later), 



(E(y)) 



dy [fJiT ■ coshy] • p(y) = U = — ■ M. (2) 

y 1 V 
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1 Accounting for the fact that multiplicity distribution of 
observed secondaries are not Poissonian 0. 



Whereas in 0] one uses Tsallis entropy instead Shannon 
ones and defines constraints in slightly different way, 
the EGE approach |10j simply adds one more constraint 
to J5J in the form of a priori known fluctuations of mean 
energy of given secondary given by its variance W, 



([E(y)-uf 



dy[vT-coshy-U} 2 -p(y) = W. (3) 
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In this case ^H] 

1 r 

p(y) = ■ cxp \~(3 ■ \i T cosh y - 7 • (U - fi T coshy)' 

(4) 

where Z is normalization constant and (3 = (3(U, W, N, fir), 
7 = "f(U,W,N, fix) are two Lagrange multipliers for the 
constraints J2J) and @ , respectively. In the case of no dy- 
namical fluctuations, i.e., 7 = 0, one recovers situation 
already known from (with some Wq = 1//3 2 with 

respect to which one should estimate effect of dynamical 
correlations) 2 . Rewriting eq. (@J as 



f 



p{y) = % ■ exp ( 



-(3* 



MT cosh y) ;/T =/3- 7- [2U-E(y)}, 

(5) 

one obtains expression formally resembling the usual Boltz- 
mann-Gibbs formula, but this time with energy-dependent 
inverse "temperature" [3* (which is thus no longer inten- 
sive variable). Actually such possibility was already dis- 
cussed in in the context of reservoir with finite heat 
capacity. It was argued there that if 
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(3{E) 
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where q is some constant, then the corresponding distribu- 
tion (where E{y) = /itrcoshy) takes form of the so called 
Tsallis distribution [H], 



v q {y) 



1 



1 



Z q (M,N) l " - *)P q (M,N).E {y)] ^, 

(7) 

with q given by © 3 . In our case where (3 = (3* one 
gets formally energy dependent Tsallis nonextensivity q 
parameter 



1 



7 



[{l3-2 1 U)+ 1 Ef 



(8) 



For 7 > it becomes smaller than unity and exceeds unity 
for 7 < 0. It coincides with result of JUj only if \*y(E — 
2U)/f3\ « 1 in which case q = 1 - 7//3 2 . 



3 Inelasticity distributions x(K) in EGE 

We have tried to apply EGE distribution as defined by 
eq. (@J to analyze the same multiparticle data as in 0] 
only to discover that these data do not require EGE, the 
best fit is obtained with 7 = or slightly negative (in 
which case the respective q from (JHJ exceeds unity, as has 
been found in 0]). The reason for this is obvious when 
inspecting Fig. ^ which confronts rapidity distributions 

2 In the center of mass frame y £ (—Y mi Ym) where Y m = 
In \M' (l + sjl - 4/4 /M' 2 ) /(2/xt)] and where M' = M - 
(N-2)fi T . 

3 Care must be taken here when considering signs because 
in one considers dependence of j3 on the energy of the 
reservoir, E R , and here we have energy of particle E = E tota i — 
Er. Therefore our q — 1 corresponds to 1 — q there. 




Fig. 1. (a) Examples of the most probable rapidity distribu- 
tions as given by eq. (|7J for hadronizing mass M = 100 GeV 
decaying into N secondaries of (transverse) mass tut = 0.4 
GeV each for different values of q parameter, (b) The same 
but using eq. Q with different values of fluctuations given by 
kw such that W = k\yWo where Wo are the intrinsic statisti- 
cal fluctuations present in the system when 7 = 0, i.e., where 
gaussian becomes pure exponential. 
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of the Tsallis type 10 with those obtained from EGE Q 
obtained for hadronization of some fixed mass M into dif- 
ferent number of secondaries. Results obtained using EGE 
show completely different behavior from Tsallis statistics 
approach clearly demonstrating that direct fluctuations in 
energy used in EGE (and characterized here by parame- 
ter kw such that W = k^Wo where Wo are the intrinsic 
statistical fluctuations present in the system when 7 = 0) 
are not equivalent to fluctuations described by parame- 
ter q of Tsallis' statistics 4 . This can be understood in the 
following way. In standard description of hadronization 
processes by means of IT in the Shannon form we always 
have some (mean) number of secondaries produced (N) 
with (mean) energy (E) each. Allowing for fluctuations of 
(N) results in q-statistics using Tsallis entropy for IT 4 . 
In this case the mean energy per particle fluctuates from 
event to event. Keeping now (N) fixed but introducing 
distribution of energy per particle (i.e., describing energy 
per particle by its mean and deviation from the mean) 
results in EGE 5 . Evidently single particle distributions in 
hadronization processes follow first or second scenario, not 
EGE. 

On the other hand EGE turns out to be very useful 
when applied to other characteristic of multiparticle pro- 
duction, namely to inelasticity distribution, x(K), (i.e., 
distribution of the fraction of the available energy, which 
is transformed into observed secondaries). In it was de- 
duced from data for the first time for two energies: 200 and 
900 GeV, cf . Fig. [21 (by analysing rapidity distributions of 
secondaries in fixed multiplicity bins). Its shape has been 
then fitted by gaussian and lorentzian curves but no expla- 
nation was offered for their possible origin and there was 
no argument at that time in favor of any of them. EGE 
provides arguments that most probably x{K) should be of 
gaussian shape. To show this let us again follow ^U] and 
let us suppose that the whole energy available for a given 
multiparticle production reaction, E = yfs, is divided into 
two parts: one part equal to E\ = K ■ y^s is going into sys- 
tem producing observed secondaries whereas the rest of it, 
E2 = E — Ei, is not used for this purpose and, in a sense, 
acts as a kind of "heath bath" (or environment) for the 
first one. Both systems, the one producing particles with 
energy E\ and the environment with energy E 2 can be in 
many possible states. Therefore 



Pi(#i) = 



f2 1 (E 1 )f2 2 (E 2 ) 



(9) 



4 Notice that for kw < 1 one gets 7 > (actually 7 — > +00 
for kw —> 0) whereas for kw > 1 one obtains 7 < leading to 
equivalent q calculated according to eq. © exceeding unity but 
otherwise being uncompatible with nonextensivity parameter 
used in upper panel of Fig. 

5 Notice that decreasing fluctuations in energy in comparison 
to standard ones, i.e., assuming in Fig. (lower panel) kw < 1 
results in tendency of particles to condensate in a single energy 
state with energy equal to E to tai/Ntotai- It means that EGE 
interpolates in fact between the microcanonical and canonical 
distributions |12|. 
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Fig. 2. Inelasticity distributions x(K) (normalized to unity) 
obtained in jl] from analysis of multiparticle production data 
for yfs = 200 GeV and yfs = 900 GeV fitted by gaussian, 
X (K) ~ exp [-(K- {K)) 2 /(2a 2 )] (full lines), and lorentzian, 
X(K) ~ cr/[A(K - (K)) 2 + a 2 ] (dash-dotted lines) formulas, 
respectively. The values of parameters (< K >;cr) for gaus- 
sian case are (0.52; 0.24) and (0.38; 0.20) for 200 GeV and 900 
GeV, respectively; for lorentzian case they are, respectively, 
(0.52; 0.25) and (0.39; 0.17) (see ;4 for more details). 



where fl denote the corresponding number of states. Defin- 
ing entropy in the usual way as 



S i (E i )=]nf2 i (E i ), 



1,2 



(10) 



one gets 



Pi(^i) 



e X p[S 1 (E 1 ) + S 2 {E 2 )]. (11) 



Expanding now entropy around E\ = U, keeping only 
linear and quadratic terms and assuming that (3 = ^- = 



din a 



E 1 = U 



and 7 = 



d z in n 

dE'f 



are the same for both 



parts of the system (generalization is straightforward) one 
immediately obtains gaussian-like form for energy E\ dis- 
tribution, 



Pi (Si) 



Z G 



exp 



- 7 (#!-[/)' 



(12) 



which, because E\ = yfs-K and U = yfs-{K) (where y/s is 
energy of reaction), translates in natural way to gaussian 
distribution of inelasticity, x(i^), as the most probably 
form with a 1 = From 0] one can also deduce the K 
dependence of the temperature T. Notice that in our case 
parameter 7 can be connected with T and heat capacity 
Cy because 7 = 1/(CVTq). On the other hand 0] pro- 
vides us also with parameter q for different inelasticities 
for both energies and we know that Cy can be connected 
with the nonextensivity parameter q, namely 1/Cy = q—1 
7 6 . It means therefore that there is simple relation con- 
necting the width a of the observed gaussian distribution 

6 Notice, however, that now CV and q are different from those 
in as they refer to both longitudinal and transverse degrees 
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Fig. 3. Inelasticity dependence of the temperature parame- 
ter as deduced in |lj from multiparticle data at energies 200 
and 900 GeV. Black circles correspond to the corresponding 
temperature obtained from the thermodynamical analysis per- 
formed along the EGE ideas. 



(i.e., parameter 7 of EGE), the temperature To and the 
nonextensivity parameter q describing internal behavior 
of the selected subsystem, namely 



T = y/2(q — 1)cte where oe = & • \A- (13) 

As can be seen in Fig. T) deduced in such way agrees 
very well with the (T, K) dependence deduced from anal- 
ysis of rapidity distributions in fixed multiplicity bins 0] 



4 Conclusions 



We would like to conclude with the following remarks: 

— EGE works only for the whole system, not for a sin- 
gle particle. This is going to be emitted according to 
its own distribution, in particular Boltzmann-Gibbs or 



of freedom, cf. |F[ for discussion how qi, and qr can be com- 
posed to produce total q. 



Tsallis, notwithstanding what the energy E\ is and 
how it is distributed. For a moment we cannot offer 
any convincing explanation why it is so. 

— EGE is not the same (i.e., it does not describes the 
same kind of fluctuations) as g-statistics. It means that 
even if for some limiting cases both distribution can be 
similar this is just an artifact. 

— On the other hand EGE tells us that for the sys- 
tem under consideration T = T(E) and E\ fluctuates. 
This means that for particles emitted from this system 
one should rather use Tsallis distributions reserving 
Boltzmann-Gibbs ones only to the case of T=const. 

Let us close with remark that the lorentzian curve 
shown also in Fig. [21 (and fitting data at least as well as 
the gaussian one) could be explained as a kind of a nonex- 
tensive extension of EGE by noticing that in q-statistical 
approach one gets gaussian distribution for q = 1 and 
lorentzian distribution for q = 2. We shall not pursue this 
further here. 
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